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Abstract 

Using deformed shell model based on Hartree-Fock intrinsic states with ^*^Zn as an example 
and employing two realistic effective interactions, namely jj44b that produces deformed shape and 
JUN45 that generates spherical shape, it is demonstrated that the On and 2iy nuclear transition 
matrix elements for double beta decay will reduce considerably as we go from spherical shapes to 
deformed shapes. This result is further substantiated by using ^^‘^Nd as another example. 
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I. INTRODUCTION 


Double beta decay (DBD) is an important open problem of current interest. There 
are many experimental efforts [1-3] and also theoretical studies ( see for example [4, 5] 
and references cited therein) in particular for neutrinoless double beta decay. The nuclear 
transition matrix elements (NTME) for neutrinoless double beta decay (0z//9/5) are 
found to vary considerably in different nuclear models. For example, QRPA values are 
always larger than the shell model values. Most variants of QRPA assume the nucleus 
to be spherical and hence they do not take into account the deformation of the nucleus. 
On the other hand, shell model takes into account all possible correlations, within the 
defined space, through the effective interaction. There are many studies which show that 
quadrupole deformation reduces the . Some of them are as follows, (i) energy density 
functional method including deformation and pairing fluctuations within the finite range 
density dependent Gogny force [6] where the authors show that deformation effects vary 
from 11% in ®^Se to as high as 38% in ^^®Te. (ii) The proton-neutron QRPA with a realistic 
residual interaction [7]. (iii) Nuclear shell model [8] where a quadrupole-quadrupole term is 
added to the effective interaction and the deformation effects were studied by varying the 
strength of the quadrupole term. In addition, there is also an attempt to study the effect of 
hexadecapole deformation [9]. In this present report, we study the effect of changing from 
spherical to deformed shape on and 2v half-lives, using realistic interactions within the 
deformed shell model (DSM) based on Hartree-Fock (HF) states. In the recent past, we 
have used DSM in the study DBD (both 2v and Oz/) in a number of nuclei in the A=60 — 90 
region [10-16]. 

In Section 2 we will discuss briefly the formulas for Oz/ and 2v half-lives and the DSM 
model. Sections 3 and 4 give the results obtained using spherical and deformed shapes for 
^°Zn and ^^°Nd nuclei respectively. Section 5 gives conclusions. 
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II. DOUBLE BETA DECAY HALF-LIVES AND DEFORMED SHELL MODEL 


Half-life for Oz/DBD for the 0^ ground state (gs) of a initial even-even nucleus decaying 
to the Oj gs of the hnal even-even nucleus is given by [4] 

[JfJM ^ Oj-)]-' = G»- |Af"‘'(0+)r , (1) 

\ / 

where {rriy) is the effective neutrino mass and is the phase space integral (kinematical 
factor). The NTME is essentially a sum of a Gamow-Teller like {Mqt) and Fermi like 
(Afp) two-body operators (we are neglecting other parts as in [15]). Then, 

M°'^(0+) = M°^(0+) - ^M°'^(0+) = (0+ II 0{2 : Oz/) II 0+) . (2) 

9a 

As seen from Eq. (2), Oz/DBD half-lives are generated by the two-body transition operator 
0{2 : Oz/). The qa and gy are the weak axial-vector and vector coupling constants. In 
our DSM analysis, we will use 0{2 : Oz/) as given in [15, 16] and it contains the neutrino 
potential, Jastrow factor and a scale factor three. It is important to mention that the results 
in the present paper are independent of the scale factor used in [15, 16] as we are studying 
only the ratio of for spherical and deformed shapes. Eqs. (1) and (2) are a result of 
the closure approximation and this is justihed only if the exchange of the light Majorana 
neutrino is indeed the mechanism responsible for the Oi/DBD. 

For 2z/ DBD, half-life for the 0^ 0^ decay is given by [19] 

[Ti/2r' = G2. |M2,|' . (3) 

The kinematical factor G 2 u is independent of nuclear structure and the M 2 u is NTME given 
by, 

^ ^ 'sp (Op II q~r+ II ij^) (ij) II II 0|) 

V [EN-{Er + EF)/2]/me ' 

Here, jOj), |0^) and |1^> are the initial, hnal and virtual intermediate states respectively 
and Em are the energies of the intermediate nucleus. Similarly, Ej and Ep are the ground 
state energies of the parent and daughter nuclei (see [12, 15] for further details). 

We will use DSM to calculate T‘^j 2 using the known values for G 2 u aud the NTME 
for nuclei with close to spherical shape and with deformation. In DSM, for a given nucleus, 
starting with a model space consisting of a given set of single particle (sp) orbitals and ef¬ 
fective two-body Hamiltonian, the lowest energy intrinsic states are obtained by solving the 
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Hartree-Fock (HF) single particle equation self-consistently. Excited intrinsic configurations 
are obtained by making particle-hole excitations over the lowest intrinsic state. These intrin¬ 
sic states do not have definite angular momenta and hence states of good angular momentum 
are projected from them. These good angular momentum states projected from different 
intrinsic states are not in general orthogonal to each other. Hence they are orthonormalized 
and then band mixing calculations are performed. DSM is well established to be a success¬ 
ful model for transitional nuclei (with A=60-90) when sufficiently large number of intrinsic 
states are included in the band mixing calculations. For details see [15] and references there 
in. 


III. RESULTS FOR SPHERICAL AND DEFORMED SHAPES FOR ™ZN 

In the DSM calculations for Oz/ and 2v NTME for ™Zn, used is the model space consisting 
of the orbitals ^ 3 / 2 , ^/ 5 / 2 , V 1/2 and ^ 5 ^ 9/2 with ^®Ni as the core. It is important to add that 
the same model space was used in many recent shell model and interacting boson model 
calculations [4, 17, 18] for ™Ge and ®^Se nuclei and ™Zn is in the same region. Within this 
space we have employed two effective interactions. These are the jj44b interaction given 
in [20] with spherical single particle energies for the four orbits being —9.6566, —9.2859, 
—8.2695 and —5.8944 MeV respectively and the JUN45 interaction given in [21] with spher¬ 
ical single particle energies for the four orbits being —9.8280, —8.7087, —7.8388 and —6.2617 
MeV respectively. With these, we have first carried out axially symmetric HF calculations 
for the parent and the daughter nuclei ™Zn and ^°Ge (for 2z/DBD, also for the intermediate 
odd-odd nucleus ™Ga) for each of the two effective interactions. 

The lowest energy HF solutions obtained using jj44b and JUN45 interactions for ™Zn are 
shown in Fig. 1 and similarly for ™Ge in Fig. 2. It is clearly seen from the HF spectra in Figs. 
1 and 2 that jj446 generates deformed solutions for the ground state of ™Zn and ™Ge. Note 
that the mass quadruple moment given in the figure gives a measure of deformation. With 
jj44b interaction, the DSM calculated energy spectra for these two nuclei are compressed 
compared to the experimental spectrum. For example, for ™Zn the calculated E{2^) is 0.355 
MeV compared to the experimental value 0.885 MeV and similarly for ™Ge they are 0.64 
MeV and 1.04 MeV respectively [16]. But, at the same time the calculated B{E2 : 2+ —)■ Oj*") 
(14 W.u. and 18.8 W.u. for ™Zn and ™Ga respectively) are in close agreement with data 
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(16.5 W.u. and 20.9 W.u. respectively). These results are reported in [16]. On the other 
hand, as seen from Figs. 1 and 2, the JUN45 interaction generates nearly spherical solutions 
for these two nuclei. For this interaction, the spectrum is in close agreement with data but 
the B{E2 : 2^ —)■ 0^) is much smaller than data value; see Fig. 3. Both the calculated energy 
spectra and B{E2) values show that jj44b generates deformed shape and JUN45 generates 
spherical shape. Therefore, calculations of NTME for Oz/ and 2u using these interactions 
will give information about shape effects on NTME. 

A. and for ™Zn 

In the DSM calculations of Ou NTME using jj44b interaction, we have considered 30 
lowest intrinsic states with K = 0^ for ™Zn, 26 lowest intrinsic conhgurations with K = 
for the daughter nucleus ™Ge (up to 6 MeV excitation in both nuclei) by making particle- 
hole excitations over the lowest HF intrinsic states generated for these nuclei. We project 
out O’*" states from each of these intrinsic states and then perform band mixing calculations 
as discussed above. Similarly for JUN45 we have considered 81 lowest intrinsic states with 
K = for ^°Zn, 68 lowest intrinsic conhgurations with K = 0^ for the daughter nucleus 
^°Ge (up to 5 MeV excitation in both nuclei) by making particle-hole excitations over the 
lowest HF intrinsic states generated for these nuclei. With jj44b interaction, the NTME 
comes out to be 2.13. However, for JUN45 interaction the value is 4.56. Thus we 
see a reduction by a factor 2 in as we go from spherical to deformed shape. Thus 

DSM is consistent with other theoretical calculations mentioned in Section 1. For further 
conhrmation of this result we will turn 2v NTME. 

For 2v DBD, hrstly as discussed in [16], [Ejv — {Ei + Ep)/2] = [1.1537 -|- Ei+] MeV and 
DSM is used to calculate Ei+ for the intermediate nucleus ™Ga as well as the numerator in 
the sum in Eq. (4). In the DSM calculations using jj44b and JUN45 interactions, we have 
considered the same intrinsic states for ™Zn and ^'^Ge as in the calculation of Oz/ NTME 
given above. In addition, for the intermediate nucleus ™Ga, 65 intrinsic states with K = 
or A = O’*" are used for jj44b interaction and 114 intrinsic states with JUN45 interaction. 
These are generated by making particle-hole excitations over the lowest HF intrinsic state 
generated for this nucleus. We project out 1^ states from each of these intrinsic states and 
then perform a band mixing calculation as discussed above. Taking the phase space factor 
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0.32 X 10“^^ yr“^ given in [19], the DSM value for is 3.39 x 10^^ yr for jj44b interaction 
and 4.94 x 10^^ yr for JUN45 interaction. Thus, for 2z/ DBD also the NTME (proportional 
to inverse of the half-life) is reduced by a factor of ~ 2.6 as we go from spherical (JUN45) 
to deformed (jj44b) shape. 


IV. RESULTS FOR SPHERICAL AND DEFORMED SHAPES FOR 

We then considered the double beta decay of ^^°Nd as another example. In the calculation 
of ^^°Nd and the daughter nucleus ^^°Sm, we used modihed Kuo-Herling interaction [23] with 
a Z = 50, N = 82 core. The protons are in the ^si/ 2 , ^^ 3 / 2 , ^^ 5 / 2 , ^( 77/2 and ^/iii /2 orbits with 
spherical single particle energies -7.92783, -7.95519, -8.67553, -9.59581 and -6.83792 MeV, 
respectively. The neutrons are in the ^pi/ 2 , ^^ 3 / 2 , ^/ 5 / 2 , ^/ 7 / 2 , ^^ 9/2 and ^* 13/2 orbits with 
spherical single particle energies for these orbits -1.050, -1.625, -1.964, -2.380, -0.895 and 
-0.261 MeV, respectively. The model space used here is same as in the recent interacting 
boson model calculations [4, 18]. As described before, we obtain the lowest energy HF 
conhguration by solving the axially symmetric HF equation self-considently. The lowest 
energy HF solutions for ^^°Nd and ^^°Sm obtained with this interaction are shown in hgures 
4 and 5. As seen from the quadrupole moments, both the nuclei have large deformations. 
To generate HF intrinsic states with nearly spherical shape, we artihcially decreased the 
effective interation to 1/lOth of its value and perform axially symmetric HF calculation as 
described above. These nearly spherical HF single particle states are shown in hgures 4 
and 5. As expected, these two solutions have small quadrupole moment and hence small 
deformation. We then proceed to calculate Oz/ NTMF with the deformed and spherical 
solutions, separately. For the deformed solutions, we considered 38 intrinsic states with 
A = O’*" for ^^°Nd. Similarly we consider 17 intrinsic states with A = O"*" for the daughter 
nucleus ^^°Sm. Now, the NTMF comes out to be 1.35. Then in the second case where 
we reduced the effective interaction matrix elements by a factor of 10 , we considered 62 
intrinsic states both for the parent and the daughter nuclei. The NTMF comes out to 
be 3.06. Thus in this case also the Oz/ nuclear matrix element increases by a factor of 2.33 
as we go from deformed to spherical shape. 
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V. CONCLUSIONS 


We have demonstrated using DSM model and ™Zn as the example and employing two 
realistic effective interactions, namely jj44b that produces deformed shape and JUN45 that 
generates spherical shape, that the Oz/ and 2v NTME for double beta decay will reduce 
considerably as we go from spherical shapes to deformed shapes. This result is further 
conhrmed by a calculation of Oz/ NTME for ^^°Nd. Our results reported in this paper are 
an important addition, since the calculations are performed using two realistic interactions 
proposed for the same model space and DSM is transperent in determining shapes, to the 
existing analysis of deformation effects on NTME for DBD. 

We would like to mention that since we have not included all the spin-orbit partners in the 
basis space, the Ikeda sum rule is not satished in our calculations. However, as mentioned 
before, many recent shell model and interacting boson model calculations [4, 17, 18] used 
the same model space. Violation of Ikeda sum rule may not be a serious dehciency in our 
calculations since we are only studying the change in double beta decay matrix elements 
with the change of nuclear shape. Calculations using much larger set of single particle orbits 
(these also need new effective interactions) are for future. 
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FIG. 1: HF single particle spectra for corresponding to lowest prolate configurations generated 
by jj44b and JUN45 interactions. In the figures circles represent protons and crosses represent 
neutrons. The Hartree-Fock energy (E) in MeV, mass quadrupole moment (Q) in units of the 
square of the oscillator length parameter and the total K quantum number of the lowest intrinsic 
states are given in the figure. Each occupied single particle orbital is two fold degenerate because 
of time reversal symmetry. 
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FIG. 2: Same as Fig. 1 but for ’^^Ge. 
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FIG. 3; The calculated yrast levels, using JUN45 interaction, for ^'^Zn and ^*^Ge compared with 
experiment. The experimental data are from ref [22]. The quantities near the arrows represent 
B(E2) values in W.u. and the effective charges used are same as in [16]. 
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FIG. 4: Same as in Fig. 1 but for with Kuo-Herling interaction (see text for details). In the 

spherical solution, the levels are degenerate as expected. 
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FIG. 5: Same as in Fig. 1 but for with Kuo-Herling interaction (see text for details). In the 

spherical solution, the levels are degenerate as expected. 
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